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Measurements of Bose-Einstein correlations of 7r~7r~ pairs in central Pb+Pb collisions were per- 
formed with the NA49 detector at the CERN SPS for beam energies of 20A, 30A, 40A, 80A, and 
158j4 GeV. Correlation functions were measured in the longitudinally co-moving "out-side-long" 
reference frame as a function of rapidity and transverse momentum in the forward hemisphere of 
the reaction. Radius and correlation strength parameters were obtained from fits of a Gaussian 
parametrization. The results show a decrease of the radius parameters with increasing transverse 
momentum characteristic of strong radial flow in the pion source. No striking dependence on pion- 
pair rapidity or beam energy is observed. Static and dynamic properties of the pion source are 
obtained from simultaneous fits with a blast-wave model to radius parameters and midrapidity 
transverse momentum spectra. Predictions of hydrodynamic and microscopic models of Pb+Pb 
collisions are discussed. 

(Septemer 26, 2007) 



PACS numbers: 25.75.-q,25.75.Gz 
I. INTRODUCTION 

In central collisions of ultra-relativistic heavy ions an 
extended volume of high energy density is created. Sev- 



eral observations at top SPS [1] and RHIC [2-5] energies 
indicate that the initial energy density is large enough 
to overcome QCD confinement and to form a phase of 
deconfined quarks and giuons during the early stages of 
the collisions. 
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This strongly compressed interacting matter expands 
into the surrounding vacuum, causing the temperature 
and energy density to drop to the critical values and forc- 
ing the system to form hadrons at the phase boundary. 
During the subsequent stage the hadrons may continue to 
scatter inelastically and elastically until all interactions 
cease. 

The excitation functions of several observables, e.g. 
the energy dependence of the pion yield per wounded nu- 
cleon and the strange particle per pion ratio indicate that 
deconfinement indeed appears already in central Pb+Pb 
collisions at lower SPS energies [6]. The transient ex- 
istence of a deconfined phase in the early stage of the 
reaction may be still reflected in such properties as e.g. 
the lifetime, final size, and collective flow of the freeze-out 
stage of the reaction [7, 8] . 

The study of Bose-Einstein (BE) correlations has been 
shown to contribute unique information on the static 
and dynamic properties of the system at thermal freeze- 
out [9, 10]. Traditionally a multi-dimensional Gaussian 
function, parametrized by radius parameters, is fitted 
to the BE correlation function. Quantitative parame- 
ters of the pion source can be derived from a simulta- 
neous analysis of single-particle transverse momentum 
spectra and two-particle BE radius parameters using hy- 
drodynamically inspired models describing the evolution 
of the system. Detailed measurements of the rapidity and 
transverse momentum dependence of BE correlations for 
7r~7r~ pairs in the full SPS energy range from 20^4 to 
IbSA GeV beam energy and their interpretation in terms 
of a rapidly expanding fireball are the subject of this pa- 
per. An alternative analysis scheme of BE correlations 
has been developed [11, 12], which tries to deduce the 
properties of the pion source directly from the BE cor- 
relations by an integral inversion technique. An analysis 
using this method has been applied to RHIC data [13] 
and is currently in progress for the SPS data of NA49. 

An earlier analysis of 7r~7r~ BE correlations at 158^ 
GeV beam energy with the NA49 spectrometer [14] pi- 
oneered the simultaneous analysis of BE correlations 
and inclusive transverse momentum spectra. The study 
found a mean velocity of the collective transverse expan- 
sion of < /? >= 0.55 ± 0.12 and a freeze-out temperature 
of T = 120 ± 12 MeV. Furthermore a strong longitudinal 
expansion and a finite duration of particle emission of 3 
to 4 fm/c were observed. The results presented in this pa- 
per basically agree with these findings. One of the major 
differences of the two analyses lies in the parametrization 
of the correlation function: in the present work the mo- 
mentum difference is decomposed into Cartesian compo- 
nents Qside, Qout, and qiong as suggested by Podgoretsky 
[15], Pratt [7], and Bertsch [16]; in the previous publica- 
tion the parametrization developed by Yano-Koonin [17] 
and Podgoretsky [15] was applied. In addition, in this 
article also measurements at lower SPS energies will be 
presented. 

Results of an analysis of K+K"*" and K~K~ BE corre- 
lations at 158^ GeV beam energy from NA49 were re- 



ported in [18]. These are less differential due to the much 
smaller number of kaon pairs. However, taking account 
of the radial flow and comparing the extracted source 
parameters at equivalent values of the transverse masses, 
consistency is found between results from pion and kaon 
analyses. 

The NA49 detector allows to study also correlations of 
heavier particles. Two-proton correlations at 158^ GeV 
beam energy were analysed in [19]. Here the contribu- 
tions from strong interactions. Coulomb repulsion and 
Fermi-Dirac quantum statistics are of similar magnitude 
and make it impossible to extract source parameters with 
the help of a simple parametrization. The final result of 
that study, i.e. the extracted averaged source radius, is 
compatible with that obtained from 7r~7r" - and K=^K=^ 
BE correlations. 

Numerous other tttt BE correlation analyses near 
midrapidity were reported in the literature. Measure- 
ments at the AGS can be found in [20] . Results from the 
SPS were pubHshed by experiments NA44 [21], WA97 
[22], WA98 [23], and NA45 [24]. Measurements at RHIC 
have been reported by STAR [25, 26], PHENIX [27] and 
PHOBOS [28]. Results show good consistency in general 
and exhibit remarkably little change over this very large 
collision energy range. At all energies, hydrodynamic 
models were found to be able to consistently describe 
many aspects of the reactions, including the transverse 
momentum dependence of BE correlations [29] . The ex- 
ception remains the over-prediction of the experimen- 
tally determined pion source radii by several fermi. This 
could be due to a too simplistic treatment of the freeze- 
out stage [30] and/or a neglect of possible refractive ef- 
fects in the fireball within the traditional interpretation 
of BE correlations [31 33]. 

This paper is organized as follows. In section 2 the 
NA49 experimental setup is described followed by a brief 
introduction to BE correlations in section 3. Selection 
criteria for events, tracks, and pairs are presented in sec- 
tion 4. Corrections, systematic uncertainties, and the 
fitting procedure are described in section 5. Results are 
presented in section 6. The paper is closed by a discus- 
sion in section 7 and conclusions in section 8. 



II. EXPERIMENTAL SETUP 

The layout of the essential components of the NA49 
detector is displayed in Fig. 1. 

Trajectories of incident beam particles are measured 
individually by beam position detectors BPD 1-3. Col- 
lisions occur in a thin lead foil T of 224 mg/cm^ thick- 
ness. The main tracking devices for charged particles 
are four large Time Projection Chambers (TPC), la- 
beled VTPCl, VTPC2, MTPCL, and MTPCR in Fig. 1. 
Two of them, VTPCl and VTPC2, are mounted in pre- 
cisely mapped magnetic fields provided by two super- 
conducting magnets with total bending power of up to 
9 Tm. The TPCs provide large acceptance, efficient 
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FIG. 1: (Color online) The NA49 experiment; beam position detectors BPD, target foil T, time of flight detectors TOF, time 
projection chambers TPC, and calorimeters CAL. Details can be found in [34]. 



reconstruction and precise tracking for charged parti- 
cles. The momentum of a charged particle can be de- 
termined from its characteristic deflection in the mag- 
netic field; with the NA49 setup a resolution of Ap/p^ w 
(0.3 — 7)-10~'* (GeV/c)~^ is achieved. A veto calorimeter 
VCAL behind a suitably adjusted collimator COLL mea- 
sures the energy of the projectile spectators. The event 
trigger consists of a valid beam track signal in coinci- 
dence with an energy deposition in the veto calorimeter 
below a preset threshold. Further details of the detector 
and trigger system can be found in [34]. The detection 
efficiency depends on particle momentum and varies from 
about 90% at midrapidity to near 100% in the forward 
direction. 

The strength of the magnetic field was scaled according 



204 GeV 1 58A GeV 




FIG. 2: (Color online) Raw yields of particle pairs at 20A and 
158A GeV beam energy as a function of and kt- 

to beam energy, i.e. at 20 A GeV it was set to «l/8 of 
the maximum setting at 158^ GeV. In this way, the ge- 



ometric acceptance in the center-of-mass system is kept 
similar at different beam energies. Fig. 2 displays the 
raw yields of particle pairs at 20^1 GeV and at 158^ GeV 
beam energy. 

The forward hemisphere in the center of mass of the 
reaction is well covered allowing differential studies in 
pion-pair rapidity 

Y.. - bn flf^'^'^^'^- (1) 

2 El + E2 - Pz.l - Pz.2 

calculated in the reaction cm. system and averaged 
transverse momentum 

fc« = ^|pt,i +Pt.2|, (2) 

where Ei and pi represent energies and momenta of the 
pions. Track topology and momentum resolution change 
with transverse momentum, rapidity, and different set- 
tings of the magnetic field. The impact of the finite 
momentum resolution on the correlation function is dis- 
cussed in detail in section VA. 

Along with their trajectory, the NA49 TPCs measure 
the specific energy loss of charged particles traversing the 
detector gas. Combined with the reconstructed momen- 
tum this information allows to identify the particle type. 
This selection is not applied in the present analysis since 
the procedure works only in a limited momentum range. 
The impact of the missing particle identification on the 
correlation function is discussed in section VB. 

The energy deposition in the veto calorimeter corre- 
sponds to the total energy of the projectile spectators 
and is closely related to the centrality of the collision. 
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At lower beam energies the threshold of the online trig- 
ger was set to accept the 7.2% most central events. The 
dataset at 158^ GeV was taken with a 10% online trigger 
setting, but to ease the comparison the event selection 
was restricted offline to the 7.2% most central events. 
The average number of wounded nuclcons for the 7.2% 
most central collisions was estimated by Glauber model 
calculations and amounts to (Nw) = 349 with a disper- 
sion of 28. 

In the course of the CERN SPS energy scan pro- 
gramme the NA49 experiment used ^°^Pb beams with 
energies of 20A, iOA, 80A, and 158^1 GcV. A sum- 
mary of the central Pb-|-Pb data sets used in this analysis 
is given in Table I. 

TABLE 1: Data sets used in this analysis. 



Beam 
Energy 
(GeV) 


(GeV) 


Year of 
data 
taking 


Number of 
analyzed 
events 


Mag. 
field 
Polarity 


20 


6.3 


2002 


360k 


+ 


30 


7.6 


2002 


420k 


+ 


40 


8.7 


2001 


217k 




40 


8.7 


2000 


360k 


+ 


80 


12.5 


2001 


296k 




158 


17.3 


1996 


386k 


+ 


158 


17.3 


2000 


502k 





III. BOSE-EINSTEIN CORRELATION 
FUNCTION 

BE correlations are a consequence of the requirement 
that the quantum mechanical wave function of two iden- 
tical bosons must be symmetric under particle exchange. 
This results in an increased probability Pi, 2 of emitting 
pion pairs at small momentum difference q = pi — P2 
relative to the product of probabilities Pi ■ P2 to emit 
single pious of momentum pi. The enhancement is usu- 
ally studied in terms of the correlation function: 



C{q) 



Pi,: 



Pi ■ P2 B{q) ■ 



(3) 



Experimentally the correlation function is constructed 

as the ratio of the pair momentum difference distribution 
S{q) and a mixed-event pair distribution B{q). This ref- 
erence distribution is constructed by forming pairs of par- 
ticles from different events. Except for genuine particle 
correlations, B{q) thus includes the features of particle 
detection with the NA49 apparatus as well as the char- 
acteristics of inclusive pion production from the source. 
The track selection criteria, described in the next sec- 
tion, are applied in the construction of the reference in 
the same way as in the construction of the signal distri- 
bution. 

The distribution of spatial distances of the emission 
points of the two pions is related to the distribution of 



their momentum difference. In fact, the spatial exten- 
sion of the pion emitting source is reflected in the width 
of the BE enhancement in the q-dependent correlation 
function. In this paper, the inverse of the width of the 
enhancement is labelled the BE radius. Elsewhere it is 
sometimes also referred to as HBT radius, in reference to 
the pioneering work of Hanbury Brown and Twiss [35] . 

Since momenta and emission points of particles from 
the source arc also correlated due to radial and longi- 
tudinal collective flow, the measured radii correspond to 
lengths of homogeneity [36] rather than to the actual ex- 
tension of the source. This feature allows to extract dy- 
namic parameters of the source, if the correlation func- 
tion is analyzed as a function of the kinematic variables 
l^,r and kf. 

The momentum difference is decomposed into 3 inde- 
pendent components as suggested by Podgoretsky [15], 
Pratt [7], and Bertsch [8]. The first, qiong = Pz,i — Pz,2, 
is the momentum difference along the direction of the 
beam, measured in the longitudinal co-moving system, 
i.e. in the frame of reference where Pz,i = ~Pz,2- The 
two other components are defined in the plane transverse 
to the beam, with qout parallel to the pair transverse 
momentum vector kt = |(Pt,i + Pt,2) and qside perpen- 
dicular to qout and qiong- 

The two-particle BE correlation function Cbe is in 
general well approximated by a three dimensional Gauss 
function 



CBB(q) = 1 -I- A • exp{ 



■2 T>2 
'Jout^out 



I? r2 

^long long 
'^QoutQlongRoutlong)- 



(4) 



A fit to the measured correlation function yields results 
for the variances Rout, Rside, Riong, and Routiong, which 
can then be compared to results from model calculations. 
For several reasons the measured correlation functions 
have amplitudes smaller than the theoretically expected 
value of two for vanishing momentum difference. The 
known sources of the suppression, such as contamination 
by misidentified particles or by weak decay products, are 
taken into account by introducing a purity factor p, de- 
scribed in section VB. To allow for additional sources of 
suppression a coherence parameter A is included in the 
fitting procedure. 



IV. EVENT, TRACK, AND PAIR SELECTION 

A. Event cuts 

Events were accepted if they had a successfully recon- 
structed beam particle and interaction vertex and an en- 
ergy deposition in the veto calorimeter within the ap- 
propriate range. The position of the main vertex was 
determined from the set of reconstructed tracks. Its po- 
sition had to be close to that of the beam track in the 
transverse plane and close to the known position of the 
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target foil in the fongitudinal direction. The deviations 
of the main vertex coordinates from the expected posi- 
tions show Gaussian distributions in x (bending plane), 
y (vertical), and z (longitudinal) directions. Events are 
accepted if their reconstructed interaction points are not 
farther away than three standard deviations from the re- 
spective mean values. 



B. Track cuts 

In this analysis, only negatively charged particles sat- 
isfying a set of quality criteria on measured track length 
and distance of closest approach to the event vertex were 
used. 

Usually tracks with less than a minimum number of 
measured points are excluded from physics analyses of 
NA49 data. However, it was found that in the VTPCl 
detector tracking quality depends on the occupancy of 
the TPCs which may vary even for fixed multiplicity in an 
event due to spiraling electrons generating large numbers 
of clusters in a TPC. These kinds of fluctuations of the 
performance of the detector cause a subtle dependence 
of the correlation function on the minimum number of 
required measured points. In order to avoid this problem 
the track quality was ensured by requiring a minimum 
number of 20 potential points simultaneously in one of 
the VTPCs and one of the MTPCs out of a possible max- 
imum of 72 and 90 points respectively. The number of 
potential points is derived from a detailed simulation and 
reconstruction of particle tracks in the detector. It corre- 
sponds to the number of charge clusters detectable in the 
TPC under ideal conditions. Varying this cut over a rea- 
sonable interval had no influence on the finally extracted 
correlations. 

To minimize contributions from products of weak de- 
cays and from secondary interactions with the detector 
material only those tracks were accepted that appeared 
to originate from the main vertex. Extrapolations of 
tracks to the nominal position of the target were required 
to deviate not more than 5.0 cm in x-direction and less 
than 3.0 cm in y-dircction from the reconstructed posi- 
tion of the main vertex. 



C. Pair cuts 

1. Two track resolution cut 

Since this study focuses on particle correlations at low 
momentum difference the two-track resolution capabili- 
ties of the detectors were investigated carefully. Tracks of 
particles with small momentum difference travel closely 
together in the TPCs. If two tracks approach each other 
the clusters may start to overlap. In this case their charge 
deposition distributions are not always resolved properly 
and the ionization points are not sufficiently well recon- 
structed. In extreme cases measured points might be 



assigned to the wrong track changing the measured val- 
ues of momentum differences significantly. Particle pairs 
possibly affected by the limited two-track resolution are 
therefore excluded from the analysis. 

In preliminary analyses of NA49 data a simple anti- 
merging cut was applied, based on an "averaged two- 
track distance" . The distance of two tracks in the trans- 
verse plane was calculated at two or more fixed planes, 
e.g. at entrance and exit of a TPC, and the average of 
these values was defined as the two-track distance. Only 
pairs were accepted, which had two-track distances larger 
than a certain minimum value. A detailed study revealed, 
that this kind of cut is not able to remove all corrupted 
pairs, resulting in a drop below unity in the projection of 
the Qout component of the correlation function visible on 
the left hand side of Fig. 3. Therefore a more involved 
procedure was developed which will be described next. 




FIG. 3; (Color online) Projections of the correlation function 
on the qout component at 158^4 GeV beam energy at midra- 
pidity in the transverse momentum bin 0.3 < kt < 0.4 GeV/c. 
Left hand side: requiring a distance of closest approach DCA 
of 1.6 cm at 10 consecutive padrows, at the right hand side: 
a DCA of 2.2 cm was required at 70 padrows. 



First each particle trajectory is extrapolated through 
the TPCs and the transverse positions at the nomi- 
nal longitudinal positions of the individual padrows are 
stored. Thus, for every pair the distance at each padrow 
is readily calculable. The trajectories are required to be 
separated by more than As = 2.2 cm at more than 50 
consecutive padrows, starting at the end of each TPC 
(see Fig. 4). If one of the tracks of a pair has less than 
50 potential points, the criterion must be fulfilled at all 
padrows where both tracks have a potential point. Start- 
ing from VTPCl, this criterion is tested for each TPC 
in which both tracks have at least 20 potential points. If 
the requirement is not fulfilled in all TPCs, the pair is 
rejected. The impact of this cut on the finally extracted 
Rout parameter is shown in Fig. 5. Varying this selection 
criterion can change the finally extracted parameters by 
up to 5%. 
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^rows 




FIG. 4: (Color online) Schematic description of the anti- 
merging cut: two tracks have to be separated by more than 
As = 2.2 cm at more than Urows ~ 50 consecutive padrows. 




As (cm) 

FIG. 5: (Color online) Dependence of the parameter Rout 
on the minimum separation As and on the required number 
of rows, at midrapidity at 0.3 < kt < 0.4 GeV/c for the 
158A GeV beam energy sample. 



2. Split track cut 

If not all charge clusters induced by a particle travers- 
ing a TPC are assigned to the track during the recon- 
struction of the track, it is possible that the remaining 
clusters are reconstructed as a separate track. In this 
case of "split tracks" (see Fig. 6) additional pairs of low 
momentum difference are artificially created vifhen con- 
structing the signal distribution. They alter the correla- 
tion function at low momentum difference significantly. 

The procedure to identify and remove pairs composed 
of the two segments of a split track is based on the follow- 
ing considerations: since both segments originate from 
the same particle they have approximately the same num- 
ber of potential points Upot- The sum of the measured 
points Umeasi}) of both track segments is equal to or 
smaller than the number of the potential points. There- 




FIG. 6: (Color online) Schematic description of a split track: 
a single track (dashed line) is reconstructed as two different 
tracks (solid lines). 

fore, a corrupted pair most probably has: 

+ 'FrT-^-^- 

npot(2) npot{2j 

Due to possible fluctuations in the determination of the 
potential and measured number of points the cut value 
is increased from 1.0 to 1.1, to ensure the removal of all 
split-track pairs. The criterion is checked for each TPC 
separately. 



D. Examples of correlation functions 

Correlation functions are constructed according to 
Eq. 3 as a three dimensional histogram with 10 MeV/c 
wide bins. The number of mixed events is chosen such 
that the background contains about eight times more en- 
tries than the signal, making the contribution from the 
background to the statistical error in the correlation func- 
tion negligible. Due to the large number of real pairs, the 
relative statistical error attributed to numerical values of 
parameters is small. The whole sample of particle pairs is 
subdivided into kinematic bins for a differential analysis. 
A bin width of 0.5 units in pair rapidity y^-jr is chosen, 
resulting in up to five bins in this direction, covering the 
range from central rapidity up to beam rapidity at the 
highest beam energy. In the transverse direction the fcj- 
value of the pair is used for the binning; five bins are de- 
fined with the boundary values 0.0, 0.1, 0.2, 0.3, 0.4, and 
0.6 GeV/c. With decreasing beam energy and increasing 
transverse momentum and pair rapidity the number of 
pair entries in the signal distribution drops rapidly pre- 
venting the fitting procedure from converging; hence fit 
results are not available in every kinematic bin at each 
beam energy. 

Fig. 7 shows expamles of projections of the measured 
correlation functions on the momentum difference com- 
ponents qside, Qout, qiong, with the integration range - 
50 MeV/c in both projected components. The examples 
were obtained for midrapidity, at 20 A, 40 A, and 158^ 
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FIG. 7: (Color online) Projections of correlations functions 
at midrapidity at three different beam energies at low (0.0 < 
kt < 0.1 GeV/c) and high (0.4 < kt < 0.6 GeV/c) transverse 
momentum. The projected components have been integrated 
over 50 MeV/c. The curves show the projections of the fits 
of Eq. 13 to the data. 

GeV beam energy and at low {kt < 0.1 GeV/c) and 
higher (0.4 < kt < 0.6 GeV/c) averaged pair transverse 
momenta. The solid lines represent the fits of Eq. 13 to 
the 3-dimensional histograms, projected in the same way 
as the measured correlation functions. 

Strong fluctuations of global event characteristics, e.g. 
of multiplicity or mean pair transverse momentum, can 
distort the reference distributions. In this analysis no 
indications of this kind have been observed, hence all 
events are sampled equally when constructing the mixed 
pair distributions. 



V. CORRECTIONS, FIT PROCEDURE AND 
SYSTEMATIC UNCERTAINTIES 

The measured correlation function is a convolution of 
BE correlations with several physics and detector related 
effects. 



The finite momentum resolution of the NA49 detector 
distorts the correlation function. The implications for 
the fitted parameters are discussed in section VA. The 
Bose-Einstein enhancement in the measured correlation 
function is diluted by pairs consisting of at least one par- 
ticle from a secondary interaction, from a weak decay, or 
a negatively charged particle of a different identity (K~ , 
p, e~). These effects are discussed in section VB and are 
corrected for by introducing a purity factor in the fitting 
procedure. 

The two particle correlation function contains not only 
BE correlations, but also correlations due to the Coulomb 
repulsion between two negatively charged particles. Since 
the aim is to extract correlation lengths due to the Bose- 
Einstein effect, a correction procedure is applied, which 
is described in section VC. 

Further contributions to the correlation function, e.g. 
due to the strong interaction or due to conservation laws, 
are assumed to be negligible in this analysis. Techni- 
cal aspects of the fitting procedure are discussed in sec- 
tion VD. 



A. Momentum resolution correction 

The accuracy of the momentum determination is lim- 
ited by the finite spatial resolution of the TPCs and by 
multiple Coulomb scattering of the particles with detec- 
tor material. The combined result is an uncertainty in 
the momentum determination which increases with in- 
creasing transverse momentum of the particle. For mea- 
surements at higher beam energies the decreasing contri- 
bution from multiple scattering leads to a slight improve- 
ment of the momentum determination. 

The finite momentum resolution causes a broadening 
of the correlation function and hence a reduction of the 
values of the fitted radius parameters. The size of the ef- 
fect was determined from a detailed simulation in which 
the resolution of the momentum difference determina- 
tion 6q{Sqout,Sqside,Sqio„g) for each q{qout,qside,qiong) 
bin was derived from a Monte Carlo study modeling the 
whole particle detection process and the reconstruction 
algorithm. 

The simulation allows to compare the generated mo- 
mentum difference of a particle pair q^"^" with the re- 
constructed value q^^'^. The differences in each compo- 
nent, e.g. q^St ~ louty filled into frequency distri- 
butions. The width of the frequency distribution de- 
fines the relative resolution in the respective component, 
e.g. Sqout- Since the momentum resolution depends 
on the particle momenta (track curvature in the TPCs) 
the procedure has to be performed differentially for each 
(l{qout,qside,qiong) bin in each (Y^^^kt) interval and has 
to be redone for each beam energy sample. In the worst 
case, i.e. at the lowest energy and for the qout component, 
the mean resolution {Sqout) is still ~5 MeV/c. 

The impact of the finite momentum resolution on the 
correlation function is estimated in the following way: 
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two artificial correlation functions Cideai and Cg 



are B. Missing particle identification and purity factor 



generated with tfieir ratio equal to the ratio of the "true" 
correlation function Ctme to the measured one Cmeas- 



deal 



(6) 



Multiplying the measured correlation function with the 
ratio Cideai I C smear allows to cxtract the BE parameters 
corrected for finite momentum resolution. 

Cideai = Sideal / Bideal represents the correlation func- 
tion measured by an ideal detector. It is derived from 
the BE-correlation function given by Eq. 4 assuming BE- 
parameters close to the fitted ones. Technically, the sig- 
nal distribution is generated by weighting the measured 
mixed-pair entries according to Eq. 4: Sideai = Bmeas ■ 
CsEiq^'^^) and the background distribution is simply de- 
fined as Bideal — Bmeas yielding Cideai ~ Sideal /B„ieas ■ 

In the construction of the correlation function C smear , 
the weight is calculated in the same way, but the posi- 
tion qsmear, whcrc thc entry is filled, is obtained from 
a smearing procedure reflecting the uncertainty in the 
momentum determination. This procedure results in 

C smear {qsrnear) = Sidealiq'^"'^) / Bmeasiq'^""') ■ The Smear- 
ing algorithm takes as input the momentum difference 
qgen ^t^q pg^jj. g^jjjj varies it slightly, using a random 
generator based on a Gaussian probability distribution 
with variance (5g(q, kt, i^Tnr) and mean equal to q^™ . 

The effect of momentum smearing on the extraction 
of the parameters is demonstrated in Fig. 8. Radius pa- 
rameters fitted with (triangles up) and without (triangles 
down) momentum smearing are compared. Clearly there 
is little impact on the fit parameters for the NA49 exper- 
iment due to its excellent momentum resolution. Hence 
a momentum resolution correction is not applied to the 
correlation functions used in the final fit. 
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FIG. 8: (Color online) Fitted radii with (triangles up) and 
without (triangles down) momentum semearing at midrapid- 
ity at 20A (upper) and 158A GeV (lower row) respectively. 



No particle type selection was applied since the particle 
identification capabilities of the NA49 TPCs allow a suf- 
ficiently clean particle by particle identification only for 
laboratory momenta between roughly 5 and 50 GeV/c. 
This restriction prevents particle identification especially 
at lower transverse momenta at midrapidity. In other 
NA49 analyses this deficiency is compensated by using 
additional information from Time of Flight detectors. 
However, because of their limited acceptance this infor- 
mation is not employed in this analysis. Instead the 
correlation functions are determined for all negatively 
charged particle pairs h^h~ and corrected for contam- 
inations. 

The fraction of pions in the sample of accepted par- 
ticles depends on the energy of the collision and on the 
kinematic region. Fig. 9 shows the rapidity distributions 
of pions and of the less abundant negatively charged 
kaons and antiprotons at different beam energies. The 
pions represent about 90% of the total number of nega- 
tively charged particles. 
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FIG. 9: (Color online) Rapidity distributions of negative pions 
and kaons and the antiproton yield at midrapidity at SPS 
energies. Open symbols indicate reflected values [6, 40, 41]. 
Data at 20A and 30A GeV are preliminary results. 

The dilution of the particle sample with non-pions re- 
duces the correlation function. Since correlations be- 
tween non-pions and any other particle, e.g. K~K~ or 
K^TT^, can be neglected, the missing particle identifica- 
tion leads only to a reduction of the fraction of correlated 
pairs, which is taken into account by introducing the pu- 
rity factor p: 



C„ 



^P-Ce, 



+ {1-P). 



(7) 



Cmeas corrcsponds to the measured correlation func- 
tion and Ccorr to the correlation function exclusively con- 
taining correlated pairs in the signal distribution. The 
purity p is defined as the ratio of the number of primary 
7r~ pairs to all negatively charged particle pairs in each 
bin of Ytttt and kt . 
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FIG. 10: (Color online) kt dependence of fit parameters at 
midrapidity at 158^ GeV with and without particle identifi- 
cation. 



FIG. 12: (Color online) kt dependence of fit parameters at 
midrapidity at 158^4 GeV assuming different values for the 
purity factor. 



This procedure has been validated by comparing ex- 
tracted parameters for h~h~ and identified 7r~7r~ pairs. 
To enable the particle identification by specific energy 
loss in the TPC the particle momenta in the laboratory 
were required to fall into the interval from 5 to 50 GeV/c. 
This selection allows a comparison only in certain kine- 
matic regions. Where the comparison is possible good 
agreement was found. An example is given in Fig. 10, 
where the kt dependence of fit parameters at midrapidity 
is shown at 158^ GeV beam energy from analyses with 
and without particle identification (PID), in the latter 
case with corrections according to Eq. 7. 
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to the total number of pion pairs is estimated from sim- 
ulations and also included in the purity factor p. Fig. 11 
shows the purity p in the analyzed y^-^ , kt bins at 158^ 
GeV beam energy. 

Fig. 12 demonstrates that the fitted radii change only 
little when instead of the purity factors shown in Fig. 11 
fixed values of 0.5 and 1.0 are used. In the following 
analyses the Yt^t^ and kt dependent purity estimates are 
used. 



C. Parametrization of the Coulomb interaction 

Identical particle correlations occur not only due to the 
quantum character of the particles. In addition, electro- 
magnetic repulsion between particles of like sign charge 
modifies the correlation function significantly at small 
momentum differences. In order to extract the correla- 
tions due to Bose-Einstein statistics, the Coulomb inter- 
action is taken into account during the fitting procedure. 
According to [42] the electromagnetic repulsion effect can 
be parametrized as a function of the modulus of the mo- 
mentum difference k* — |k*| with k* = = and 
the mean pair separation (r*), where the asterisk indi- 
cates values calculated in the c.m.s. of the particle pair. 
If the momentum difference is smaller than /c, defined as 



k = 



4(r* 



1 + 2^ l + d2 — 
a \ a 



(8) 



with the TT TT Bohr radius a — 388 fm and d2 — Stt/S, 
the Coulomb weight is given by (fc* < k): 



FIG. 11: (Color online) The purity factor p as function of 
and kt derived from simulations at 158^ GeV beam energy. 

The measured track sample is not only contaminated 
with non-pions. Also weak decays of heavier particles 
away from the vertex can generate pions which cannot 
be distinguished experimentally from primary particles. 
These decay products show no correlations with pions 
generated in the interaction volume. Their contribution 
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, (9) 



with the Gamov term Ac{ri) = 27r?7(exp(27r?7) — 1)^^ and 
T] = l/{k*a). In case k* is larger than fc, the Coulomb 
weight is given by {k* > k): 



C C oulomh 



a{r*)k*' 



(10) 
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with 



D. Fitting procedure 



d = a{r*)k l-Ac{r,) 1 + 2^ l + dz — 

L L V ^ 

' ^ (11) 

Under the assumption that both (r*) and the BE param- 
eters are fixed by the same freeze-out conditions and that 
Riong - Rside - R and = ^Rout the mean pair 
separation (r*) is approximated in [42] by 



--r: 



1,1 



2e 1 



(12) 



where e = (1 - (i?/i?Lt)')'/' 

During the fitting procedure R is approximated by i? ~ 
{Riong + Rside) /'^ and (r*) is obtained iteratively from 
Eq. 12. 

Fig. 13 shows fitted BE parameters using Coulomb 
weights based on (r*) derived from the fitted radii as 
well as results obtained with fixed values of (r*). Sys- 
tematic variations of less than 5 % are found in the radii 
for the different assumptions. 
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FIG. 13: (Color online) Fit parameters (158A GeV, midrapid- 
ity) using various values for the mean pair separation (r*) in 
the Coulomb correction procedure. Note the different scales 
on the the vertical axes. 

With the NA49 setup it is possible to measure the cor- 
relation function of oppositely charged particles. Based 
on this measurement a Coulomb correction procedure 
can be derived, not relying on any assumption about the 
mean pair separation. In this way the Coulomb interac- 
tion was taken into account in a previous NA49 publica- 
tion [14]. Since in certain kinematic regions the fraction 
of positively charged particles which do not contribute 
to the correlation function is only poorly known but cer- 
tainly relatively large, the analytical procedure is pre- 
ferred in the present analysis. 



The measured correlation function Cmeasio) 
S{q) / B{q) is fitted by the parametrization: 



Cyneasiq) ^ 71 ■ [p ■ CbE ' Ccoulomb + (1 "P)), (13) 

where n represents a normalization factor introduced to 
account for the different number of entries in the sig- 
nal and reference distributions. The normalization can 
be derived in different ways. In this analysis it is deter- 
mined in the fitting procedure where it is treated as a free 
parameter. The remaining terms of Eq. 13, the purity p, 
the Bose- Einstein correlation Cbe, and the Coulomb cor- 
relation Ccouiomb have been introduced in the preceding 
paragraphs. 

The free parameters, i.e. the radii Rside, Rout, Riong, 



Rq 



Ltlo 



the incoherence parameter A and the normali- 



sation constant n are obtained from a minimization of 



= Y^{S{q,)/B{qi)-Cmeas{qi)f/<1^l)^ (14) 



e{q^) = S{q,)/B{q,) ■ ^l/S{q,) + l/B{q,) 

using the MINUIT package [43]. The fitting routine in- 
cludes all bins i with |q| < 0.2 GeV/c. Varying the fitted 
range changed the resulting fit parameters by less than 
5%. The MINUIT package attributes an error to each 
parameter, which is equal to the inverse of the second 
derivative of the function at its minimum with re- 
spect to this parameter. Along with the numerical value 
of each parameter at the minimum of the function its 
uncertainty is always given. A quantitative estimate of 
the fit quality is not given, since it depends strongly on 
the fitted range and details of the fitting procedure. How- 
ever, comparing the projections of the fitted and mea- 
sured correlation functions in Fig. 7, it can be concluded 
that Eq. 13 describes the data very well. 



E. Systematic uncertainties 

The uncertainty of the extracted BE parameters is 
characterized not only by the error estimated during the 
fitting procedure, but also by the systematic uncertain- 
ties in the whole analysis procedure. Indeed, due to the 
large statistics of the pair samples, the systematic errors 
are expected to be larger than the statistical ones. 

In order to estimate the systematic errors, the crucial 
analysis cuts are varied over a reasonable range and the 
impact of different sets of cuts on the fit parameters is 
calculated. The following paragraphs discuss the individ- 
ual cuts which were varied and the impact on the values 
of the radii. 

A detailed comparison to the analysis presented in [24] 
revealed that the way of normalizing the correlation func- 
tion may be relevant. In this analysis the normalization 
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is treated as a free parameter in the fit procedure, but 
it can also be derived from the ratio of the number of 
entries in the signal distribution to the number of entries 
in the background distribution. The difference between 
the fit parameters determined using the two different ap- 
proaches is taken as an estimate for the systematic un- 
certainty due to the normalization method. On average, 
the systematic uncertainty due to the unknown normal- 
ization amounts to about 0.1 fm in the radii. 



The way the Coulomb interaction is taken into ac- 
count influences the final results. The size of this effect 
is estimated by assuming different values for the mean 
pair separation in Eqs. 8-11. The fit is done assuming 
(r*) = 10 fm and (r*) = 30 fm and the resulting parame- 
ters are compared to the standard procedure where (r*) is 
determined iteratively from the radius parameters. The 
maximum difference between the standard procedure and 
the two extreme cases is taken as an estimate of the sys- 
tematic uncertainty due to the treatment of the Coulomb 
interaction. Also this uncertainty may change the final 
values of the radii by the order of 0.1 fm. 



The determination of the purity factor may also influ- 
ence the final results in a systematic way. Hence the fit 
was done assuming a purity 20% larger and 20% smaller 
than the purity derived from the simulations. Again, 
the maximum difference between the standard approach 
and the extreme values is taken as an estimate or the 
uncertainty introduced by the purity determination and 
amounts in most cases to less than 0.1 fm. 



One of the most sensitive pair cuts is the anti-merging 
cut described in paragraph IV C 1. As can be seen from 
Fig. 5 the fit results change considerably with the number 
of required rows and the minimum separation. Varying 
the cut values in a reasonable range (45 < number of 
rows < 72 and 1.8 cm < minimum separation < 3.0 cm) 
and comparing the fit results to the standard cut (num- 
ber of rows = 50 and minimum separation =2.2 cm) the 
size of this systematic uncertainty is estimated. Since it 
is very time-consuming this procedure is only performed 
for the bin with the largest expected uncertainty, i.e. at 
high transverse momentum, at midrapidity and at 158^ 
GeV beam energy. The systematic error due to this un- 
certainty in the reconstruction process reaches values up 
to 0.5 fm. 



The overall systematic error assigned to the radius pa- 
rameters is defined as the maximum value of the four er- 
rors described above. For the A parameter the variation 
of the purity value was not considered, since these two 
parameters are directly correlated. In Table II both the 
error from statistics and the systematic error are given. 
The figures in section VI show statistical errors only. 



cP 



R 



side 



R 



out 



'^3 



A □ 



4 
2 



£ O D 

6- OA 



4 

2 



0.5 1 

k, (GeV/o) 



0.5 1 

k, (GeV/c) 



0.5 1 

k, (GeV/c) 



E 

^■5 

-5 



R 



outlong 



^ 1 
0.8 
0.6 
0.4 
0.2 



• NA49 
o WA98 
A NA45 
» NA44 
□ WA97 



0.5 1 

k, (GeV/c) 



0.5 1 

k, (GeV/c) 



FIG. 14: Dependence of fit parameters on kt at midrapidity 
at 158^4 GeV beam energy. Open symbols indicate measure- 
ments by other experiments [21-24]. 



VI. RESULTS 

A. Energy and kt dependence of fit parameters at 
midrapidity 

Fig. 14 shows the kt dependence of the extracted BE 
parameters at midrapidity (2.9 < Y^tt < 3.4) at the top 
SPS beam energy of 158^ GeV from this analysis (filled 
dots) along with previously published results from exper- 
iments NA44 [21], WA97 [22], WA98 [23], and NA45 [24]. 
The parameters Rside, Rout, and Riong decrease with in- 
creasing transverse momentum. In heavy-ion collisions, 
this behaviour is commonly attributed to collective ex- 
pansion, with the rate of decrease reflecting the strength 
of the expansion. In section VII, the kt dependence will 
be compared to model calculations. At vanishing trans- 
verse momentum Rout and Rside are expected to con- 
verge. Indeed the increase of Rout flattens at low kt and 
the values of Rout approach those of Rside ■ The param- 
eter Routiong remains small, independent of transverse 
momentum, as expected at midrapidity. 

The A parameter deviates from unity. Since the dilu- 
tion of the primary pion pairs by secondary pions and 
non-pions is already included via the purity factor, the 
departure of the A parameter from unity has to be at- 
tributed to physics effects. One often proposed expla- 
nation is the incomplete chaoticity of the source. The 
value of the A parameter could also be related to fea- 
tures of the pion emission, e.g. [37] discusses the impact 
of the non-zero proper volume of the pions on the cor- 
relation function. It could also be partially caused by a 
non-Gaussian shape of the correlation function due to e.g. 
long-lived resonance states, or the approximation used to 
parametrize the Coulomb interaction [38] . Details of the 
functional form of the correlation function are currently 
studied using source imaging methods [39] introduced in 
[11, 12]. 
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Fig. 14 also displays published results from other ex- 
periments [21-24] for comparison. One should note that 
the analysis procedures applied by the collaborations dif- 
fer in details, e.g. in the treatment of the Coulomb in- 
teraction. Moreover, some experiments did not use the 
mean transverse momentum variable fct, but some similar 
measure. In Fig. 14 differences between those kinematic 
variables are neglected. Absolute values as well as trans- 
verse momentum dependence of radii extracted by the 
various experiments are similar. 



lies in the treatment of the Coulomb interaction or in the 
fitting procedure. Rather it is due to differences in the 
measured correlation functions themselves. 
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FIG. 15: Dependence of fit parameters on kt at midrapidity 
for difi'erent beam energies. Full dots show results from this 
analysis, open circles indicate measurements of NA45 [24] . 
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FIG. 16: Dependence of the fit parameters on rapidity at low 
kt at different beam energies. Filled symbols correspond to 
0.0<fct < 0.1 GeV/c, open symbols to 0.1 < h < 0.2 GeV/c. 



In Fig. 15 the kt dependence of the fit parameters is 
plotted at five different beam energies in the central ra- 
pidity bin. At all energies, the transverse momentum 
dependence of the fit parameters is similar. The radii 
decrease with increasing transverse momentum, neither 
absolute values nor gradients change significantly with 
increasing beam energy. The cross term Routiong shows 
small values independent of energy and transverse mo- 
mentum. The A parameter decreases weakly with in- 
creasing beam energy. 

At 40 A and 80 A and 158 A GeV beam energy, measure- 
ments of NA45 [24] are also shown. The radii agree in 
the side and long component, but NA45 reports system- 
atically smaller values of Rout- The largest discrepancy 
is observed at 40A GeV beam energy. A detailed com- 
parison study excluded that the origin of the difference 



B. Dependence of fit parameters on rapidity 

The dependence of the fit parameters on rapidity is 
shown in Figs. 16 and 17 for the 5 beam energies of this 
analysis. At lower transverse momenta (see Fig. 16), the 
parameter Rside decreases weakly with increasing rapid- 
ity, with a larger gradient at smaller beam energies. The 
parameter Rout remains almost constant at all rapidities 
and at all energies. The rapidity dependence of the pa- 
rameter Riong shows no clear systematic trend. At small 
transverse momenta (0.0 < kt < 0.1 GeV/c) variations 
of the order of one fermi are observed. The parame- 
ter Routiong increases clearly with increasing rapidity, as 
expected for particle sources with non- vanishing longitu- 
dinal expansion [44]. The A parameter shows no distinct 
rapidity dependence, but it decreases with increasing en- 
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FIG. 17: Dependence of the fit parameters on rapidity at high 
kt at different beam energies. Filled symbols correspond to 
0.2 <kt < 0.3 GeV/c, open symbols to 0.3 < fc* < 0.4 GeV/c. 



ergy from values of about 0.8 at 20^1 GeV beam energy 
to about 0.6 at 158A GeV beam energy. 

The rapidity dependence of the fit parameters at larger 
transverse momenta is displayed in Fig. 17. Since the 
pion yield decreases with increasing transverse momen- 
tum and with increasing rapidity, the statistics in several 
kinematic bins, especially at low beam energies at larger 
rapidities, is too small to obtain a reasonable fit result. 
The parameters Rside, Rout, and at larger transverse mo- 
menta also Riong, depend only weakly on rapidity and 
on beam energy. The parameter Routiong increases with 
increasing rapidity also at larger transverse momenta. 
The A parameter shows larger fluctuations in the case 
of higher transverse momenta, but there is no obvious 
pattern visible. 



C. Numerical values 

Tables II and III list numerical values of all fit param- 
eters. The first three columns indicate beam energy, ra- 
pidity bin and mean transverse momentum respectively. 



The next columns give the numerical values of the fitted 
parameters with their statistical and systematic error in 
parentheses. For completeness the last two columns give 
the estimated purity factor and the mean pair separa- 
tion used for the Coulomb correction and calculated it- 
eratively from the fitted radius parameters. 



VII. DISCUSSION 
A. Relation between radii and source parameters 

The correlation function of identical pions measured 
in heavy-ion collisions is determined by the space-time 
structure of the system at kinetic freeze-out. The condi- 
tions of the kinetic freeze-out result from the expansion 
of the dense and strongly interacting matter created in 
the collision into the surrounding vacuum. The evolu- 
tion of this expansion is defined by the properties of the 
matter itself; in this way the equation of state of nuclear 
matter is linked to the measured correlations. 

The relation between the dynamic evolution of the sys- 
tem and the measured radii, characterizing the correla- 
tion function, cannot be given in a model-independent 
way. Nevertheless, making several general assumptions, 
a set of equations connecting key parameters of the ex- 
pansion to the magnitude of the measured radii and their 
dependence on the transverse mass have been derived 
(see e.g. [9] and references therein') . With the transverse 
mass defined as mt — {m'^ + fcj )^", where m^r indicates 
the pion mass, and (3t = kt/rrit, the approximate rela- 
tions are: 



Rlong 
Rside 



d2 
^out 



r,2 



TO • (T/mO^/^ (15) 
Rg.o/il + mt-v'j/T)^/^ (16) 

(17) 



Ar2 



Pi 



Here tq denotes the total lifetime of the system and 
T the temperature, Rgeo the transverse size of the pion 
source and r]f the strength of the transverse expansion at 
freeze-out. At represents the duration of particle emis- 
sion. Although, Eqs. 15 - 17 only hold under a number of 
assumptions, they may serve to guide the interpretation 
of the measured correlation parameters. 

The Rlong parameter is connected to the total lifetime 
of the system. In Figs. 15 - 17 the fct-dependence of 
Rlong changes only slightly with beam energy and ra- 
pidity interval. Assuming a fixed freeze-out temperature 
this observation means that the lifetime of the system is 
independent of pion pair rapidity and of incident beam 
energy. 

Also the A:(-dependence of Rside and the intercept when 
extrapolated to vanishing transverse momentum change 
little with beam energy, consequently the transverse ra- 
dius of the particle source at freeze-out and the strength 
of the transverse expansion vary only weakly with beam 
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TABLE II: Numerical results of the fitting procedure at 20A, 30A, and 40A GeV beam energy. The first three columns give 
the beam energy, the rapidity interval, and the mean transverse momentum, respectively. The next five columns show the BE 
parameters, in the last two columns the purity factor and the assumed mean pair separation applied in the Coulomb correction 
procedure are listed. 



(AGeV) 



(GeV/c) 



Rside 

(fm) 



Rout 

(fm) 



Rlong 

(fin) 



R outlong 

(fm) 



(r*) 
(fm) 



20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 



30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 
30 



40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 



0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 



0.06 
0.15 
0.25 
0.34 
0.48 



5.90 
5.46 
4.73 

4.51 
4.04 



0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 



0.06 
0.15 
0.25 
0.34 



5.70 
5.24 
4.36 
4.05 



1.0-1.5 
1.0-1.5 
1.0-1.5 



0.06 
0.15 
0.24 



5.37i 
4.95i 
4.20 



1.5-2.0 
1.5-2.0 



0.06 
0.14 



5.08 
4.73 



2.0-2.5 0.06 4.87i 



0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 



0.06 
0.15 
0.25 
0.35 
0.48 



5.88 
5.25 
4.65 
4.11 
3.66 



0.5-1.0 
0.5-1.0 
0.5-1.0 

0.5-1.0 
0.5-1.0 



0.06 
0.15 
0.25 

0.35 
0.48 



5.54 
5.001 
4.37i 
3.93i 
3.36i 



1.0-1.5 
1.0-1.5 
1.0-1.5 
1.0-1.5 



0.06 
0.15 
0.24 
0.34 



5.43i 
4.97i 
4.59i 
4.77i 



1.5-2.0 
1.5-2.0 
1.5-2.0 



0.06 
0.14 
0.24 



5.21 
4.77i 
4.48 



2.0-2.5 
2.0-2.5 



0.06 
0.14 



4.91 
4.34 



0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 



0.06 
0.15 
0.25 
0.35 
0.48 



5.64 
5.43 
4.63 
4.16 
3.73 



0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 

0.5-1.0 



0.06 
0.15 
0.25 
0.35 

0.48 



5.60 
5.30 
4.25 
3.79 

3.60 



1.0-1.5 
1.0-1.5 
1.0-1.5 
1.0-1.5 



0.06 
0.15 
0.24 
0.34 



5.41 
5.06 
4.56 
4.26 



1.5-2.0 
1.5-2.0 
1.5-2.0 
1.5-2.0 



0.06 
0.15 
0.24 
0.34 



5.22 
4.80 
4.25 
3.93 



2.0-2.5 
2.0-2.5 
2.0-2.5 



0.06 
0.14 
0.24 



4.97i 
4.48 1 
3.99 



0.09 
0.08 
0.08 
0.12 
0.14 



0.08 
0.09 
0.10 
0.16 



0.08 
0.10 
0.15 



0.11 
0.11 



0.19 



0.08 
0.07 
0.08 
0.10 
0.11 



0.06 

0.07' 
0.10 
0.13 
0.19 



0.06 
0.06 
0.12 
0.23 



0.08 
0.07 
0.13 



0.14 
0.13 



0.06 
0.06 
0.07 
0.09 
0.10 



0.06 
0.05 

0.09 
0.14 
0.20 



0.05 
0.05 
0.08 
0.16 



0.06 
0.06 
0.08 
0.16 



0.10 
0.10 
0.16 



0.16 
0.13 
0.14 
0.35 
0.79 



0.13 
0.13 
0.20 
0.59 



0.13 
0.15 
0.35 



0.15 
0.25 



0.52 



0.16 
0.13 
0.14 
0.27 
0.48 



0.13 

0.13 
0.14 
0.36 
1.00 



0.13 
0.13 
0.16 
0.61 



0.13 
0.13 
0.32 



0.24 
0.34 



0.15 
0.13 
0.14 
0.20 
0.36 



0.13 
0.13 

0.14 
0.27 
0.56 



0.13 
0.13 
0.14 
0.38 



0.13 
0.13 
0.16 
0.69 



0.13 
0.16 
0.53 



6.25 
5.95 
4.98 
4.52 
3.56 



5.95i 
6.17i 
5.I81 
3.87i 



5.85 
6.141 
5.50 



5.85 
6.02 



5.46 



6.3O1 
6.07i 
5.21 
4.72i 
4.051 



6.03 

5.95i 
5.17i 
4.38i 
3.681 



5.90 
5.94 
5.96 
3.59 



5.71 
6.30 
5.71 



5.45 
6.46 



6.32i 
6.04i 
5.37i 
4.59i 
4.02i 



5.95 
6.08 

5.72 
4.60 
5.01 



5.881 
6.29i 
5.77i 
5.OO1 



5.83 
6.25 
5.97i 
4.86 



5.72 
6.69 
5.141 



0.09 
0.08 
0.10 
0.14 
0.16 



0.08 
0.11 
0.14 
0.18 



0.09 
0.13 
0.25 



0.13 
0.15 



0.23 



0.08 
0.08 
0.10 
0.12 
0.15 



0.07 

0.09 
0.13 
0.21 
0.26 



0.07 
0.08 
0.18 
0.34 



0.09 
0.10 
0.19 



0.17 
0.21 



0.07 
0.07 
0.09 
0.12 
0.14 



0.06 
0.07 

0.14 
0.27 
0.43 



0.06 
0.07 
0.11 
0.20 



0.07 
0.08 
0.12 
0.24 



0.12 
0.16 
0.26 



0.20 
0.22 
0.35 
0.49 
0.89 



0.20 
0.25 
0.35 
0.75 



0.20 
0.22 
0.75 



0.20 
0.37 



0.66 



0.20 
0.25 
0.35 
0.36 
0.75 



0.20 

0.24 
0.35 
0.78 
1.40 



0.20 
0.22 
0.35 
1.72 



0.20 
0.27 
0.58 



0.32 
0.58 



0.20 
0.25 
0.35 
0.30 
0.54 



0.20 
0.23 

0.40 
0.73 
1.68 



0.20 
0.26 
0.35 
0.58 



0.20 
0.30 
0.35 
1.48 



0.20 
0.39 
1.30 



6.91 
5.80 
4.42 
3.53 
2.94 



6.56 
5.69 
4.19 
3.31 



7.10 
5.90 
4.32 



7.93 
6.24 



8.65 



7.47 
6.05 
4.51 
3.55 
2.85 



7.06 

5.69 
4.34 
3.45 
2.42 



7.21 
5.88 
4.66 
4.071 



7.71 
6.43 
5.071 



8.43 
6.91 



7.72 
6.41 
4.74 
3.73 
3.15 



7.19 
6.00 

4.42 
3.31 
2.73 



7.17 
6.15 
4.56 
3.69 



7.72 
6.53 
4.92 
4.18 



8.75 
7.23 
5.21 



0.11 
0.08 
0.08 
0.10 
0.14 



0.10 
0.10 
0.11 
0.14 



0.12 
0.13 
0.17" 



0.18 
0.16 



0.38 



0.10 
0.08 
0.08 
0.10 
0.09 



0.09 

0.08 
0.11 
0.13 
0.15 



0.09 
0.09 
0.14 
0.20 



0.13 
0.11 
0.16 



0.28 
0.25 



0.09 
0.07 
0.08 
0.09 
0.11 



0.08 
0.07' 
0.09 
0.15 
0.17 



0.07 
0.07 
0.09 
0.13 



0.10 
0.08 
0.11 
0.18 



0.21 
0.17 
0.24' 



0.09 
0.23 
0.24 
0.32 
0.83 



0.09 
0.23 
0.25 
0.50 



0.13 
0.23 
0.46 



0.28 
0.37 



1.21 



0.13 
0.23 
0.24 
0.30 
0.42 



0.10 

0.23 
0.24 
0.47' 
1.05 



0.13 
0.23 
0.24 
0.45 



0.26 
0.23 
0.46 



0.59 
0.71 



0.25 
0.23 
0.24 
0.22 
0.36 



0.16 
0.23 

0.24 
0.30 
0.40 



0.14 
0.23 
0.24 
0.34 



0.27 
0.23 
0.24 
0.91 



0.59 
0.36 
0.93 



1.22 
1.41 
0.00 
1.10 
0.37 



2.45 
2.51 
2.44 
1.38 



3.30 
3.76 
3.54 



4.54 
4.54 



4.86 



1.71 
1.86 
1.51 
1.18 
1.18 



2.21 

2.77 
2.52 
2.32 
1.36 



3.43 
3.61 
3.60 
1.95 



4.35 
4.86 
4.26 



4.72 
5.58 



1.88 
1.05 
0.73 
0.85 
0.00 



2.54 
2.671 
2.24 
1.32 
1.62 



3.18 
3.87 
3.43 
2.77 



4.19 
4.73 
4.24 
3.45 



5.13 
5.83 
4.16 



0.38 
0.25 
1.34 
0.28 
0.65 



0.16 
0.18 
0.17 
0.24 



0.14 
0.16 
0.23 



0.17 
0.17 



0.32 



0.24 
0.18 
0.18 
0.23 
0.20 



0.16 

0.13 
0.16 
0.19 
0.24 



0.11 
0.10 
0.18 
0.31 



0.12 
0.11 
0.18 



0.23 
0.24 



0.20 
0.28 
0.37 
0.30 
0.92 



0.12 
0.11 

0.18 
0.32 
0.30 



0.09 
0.08 
0.11 
0.19 



0.10 
0.09 
0.12 
0.21 



0.17 
0.17 
0.23 



0.08 
0.08 
0.00 
0.11 
0.74 



0.16 
0.18 
0.17 
0.19 



0.22 
0.27 
0.52 



0.32 
0.31 



0.73 



0.15 
0.13 
0.08 
0.08 
2.35 



0.18 

0.23 
0.17 
0.47 
0.76 



0.26 
0.29 
0.24 
0.94 



0.32 
0.38 
0.44 



0.40 
0.57 



0.55 
0.07 
0.06 
0.13 
0.50 



0.22 
0.25 
0.22 
0.28 
0.45 



0.26 
0.34 
0.25 
0.37 



0.35 
0.40 
0.30 
1.05 



0.50 
0.50 
1.04' 



0.79 
0.68 
0.64 
0.68 
0.68 



0.79 
0.71 
0.56 
0.53 



0.81 
0.73 
0.63 



0.83 
0.77 



0.86 



0.76 
0.67 
0.58 
0.56 
0.58 



0.771 
0.66 
0.54 
0.52 
0.38 



0.78 
0.72 
0.69 
0.99 



0.76 
0.75 
0.78 



0.76 
0.77 



0.69 
0.67 
0.57 
0.56 
0.58 



0.69 
0.68 

0.54 
0.40 
0.48 



0.73 
0.73 
0.67 
0.70 



0.75 
0.73 
0.72 
0.771 



0.77 
0.72 
0.83 



0.02 
0.01 
0.02 
0.03 
0.05 



0.01 
0.02 
0.02 
0.04 



0.02 
0.02 
0.04 



0.02 
0.03 



0.05 



0.01 
0.01 
0.02 
0.02 
0.03 



0.01 

0.01 
0.02 
0.04 
0.05 



0.01 
0.01 
0.03 
0.13 



0.02 
0.02 
0.04 



0.04 
0.04 



0.01 
0.01 
0.01 
0.02 
0.03 



0.01 
0.01 

0.02 
0.04 
0.07 



0.01 
0.01 
0.02 
0.04 



0.01 
0.01 
0.02 
0.06 



0.03 
0.03 
0.07 



0.15 
0.13 
0.09 
0.07 
0.18 



0.14 
0.11 
0.07 
0.09 



0.14 
0.12 
0.05 



0.13 
0.12 



0.11 



0.14 
0.12 
0.08 
0.06 
0.09 



0.14 
0.11 
0.06 
0.08 
0.52 



0.13 
0.12 
0.11 
0.05 



0.12 
0.11 
0.09 



0.09 
0.10 



0.13 
0.12 
0.08 
0.06 
0.08 



0.13 
0.12 
0.08 
0.05 
0.11 



0.13 
0.12 
0.11 
0.09 



0.12 
0.11 
0.09 
0.15 



0.09 
0.09 
0.07 



0.39 14.76 
0.51 15.17 
0.63 15.12 
0.67 16.43 
0.65 16.50 



0.45 14.12 

0.53 15.23 

0.62 15.14 

0.66 14.32 



0.51 14.20 
0.57 15.11 
0.65 15.73 



0.58 14.41 
0.64 15.00 



0.67 13.40 



0.41 
0.51 
0.59 
0.66 
0.64 



15.22 
15.37 
15.60 
16.77 
18.15 



0.45 14.45 

0.53 14.79 

0.60 15.23 

0.66 15.69 

0.66 16.41 



0.51 14.37 

0.55 14.86 

0.62 17.06 

0.66 14.43 



0.59 14.12 
0.63 15.51 
0.68 16.71 



0.69 13.42 
0.70 15.67 



0.43 
0.53 
0.58 
0.62 
0.62 



15.26 
15.70 
16.03 
16.54 
18.20 



0.47 14.53 
0.54 15.37 
0.57 16.30 
0.61 16.16 
0.63 21.67 



0.52 14.32 

0.55 15.56 

0.57 16.61 

0.59 17.66 



0.59 14.44 
0.62 15.56 
0.65 16.98 
0.66 17.32 



0.65 14.13 
0.68 16.27 
0.69 15.32 



energy. Closer inspection of Fig. 16 shows a system- atic decrease of the Rside parameter at low transverse 
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TABLE III: Numerical results of the fitting procedure at 80^ and 158^ GeV beam energy. The first three columns give the 
beam energy, the rapidity interval, and the mean transverse momentum, respectively. The next five columns show the BE 
parameters, in the last two columns the purity factor and the assumed mean pair separation applied in the Coulomb correction 
procedure are listed. 



Beam Y„r 
{AGeV) 



(GeV/c) 



Rside 

(fm) 



Rout 

(fm) 



Rlong 

(fm) 



R outlong 

(fm) 



(r*) 
(fm) 



80 
80 
80 

80 
80 



80 
80 
80 
80 
80 
80 
80 
80 
80 



80 
80 
80 



158 

158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 
158 



0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 



0.06 
0.15 
0.25 
0.35 
0.48 



5.73 

5.20i 
4.74i 
4.011 
3.50i 



0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 



0.06 
0.15 
0.25 
0.35 
0.48 



5.27i 
4.97i 
4.75i 
4.38i 
3.87i 



1.0-1.5 
1.0-1.5 
1.0-1.5 
1.0-1.5 

1.0-1.5 



0.06 
0.15 
0.25 
0.34 
0.47 



5.30 
5.13 
4.75 
4.24 
3.41 



1.5-2.0 
1.5-2.0 
1.5-2.0 
1.5-2.0 



0.06 
0.15 
0.24 
0.34 



5.14 
4.81 
4.14 
4.03 



2.0-2.5 
2.0-2.5 
2.0-2.5 



0.06 
0.15 
0.24 



4.66 
4.28 
4.29 



0.0-0.5 

0.0-0.5 
0.0-0.5 
0.0-0.5 
0.0-0.5 



0.06 

0.15 
0.25 
0.35 
0.48 



5.59i 

5.25i 
4.66i 
4.42i 
4.17i 



0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 
0.5-1.0 



0.06 
0.15 
0.25 
0.35 
0.48 



5.30i 
5.O61 
4.67i 
4.29 
4.08 



1.0-1.5 
1.0-1.5 

1.0-1.5 
1.0-1.5 
1.0-1.5 



0.06 
0.15 
0.25 
0.35 
0.48 



5.03 
5.11 

4.78 
4.36 
4.10 



1.5-2.0 
1.5-2.0 
1.5-2.0 
1.5-2.0 
1.5-2.0 



0.06 
0.15 
0.25 
0.34 
0.47 



5.42 
4.98 
4.68 
4.15 
3.86 



2.0-2.5 
2.0-2.5 
2.0-2.5 
2.0-2.5 



0.06 
0.15 
0.24 
0.34 



5.00 

4.79 
4.46 
4.65 



0.08 

0.07' 
0.08 
0.13 
0.23 
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momenta with increasing rapidity, suggesting a smaller 
transverse radius at larger rapidities. 

At lower kt, the parameter Rout is always larger than 
Rside, indicating a finite emission duration. In some cases 
at large transverse momenta, the radius parameter Rout 
is smaller than Rside, rendering relation 17 meaningless. 
At RHIC energies some analyses also found smaller values 
of Rout compared to Rside [25], triggering several detailed 
theoretical studies, e.g. [45]. 



Relations 15 - 17 suffer from the fact, that only two- 
particle correlation functions are considered, while other 
observables related to the collective behaviour of the 
fireball are disregarded. Single particle spectra and 
anisotropic flow measurements also carry information on 
the expansion and freeze-out temperature and should 
therefore be addressed together with the BE correlations 
in order to get a consistent and complete picture of the 
evolution of the particle-emitting source. In fact, the 
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ambiguity of transverse expansion and temperature in 
Eqs. 15 and 16 can only be overcome if additional infor- 
mation is provided, e.g. by a simultaneous fit of single 
particle spectra and the fct-dependence of the correlation 
radii, as performed for the first time in [14]. Below, the 
radii will be compared to more comprehensive model cal- 
culations, which indeed aim to describe not only the two- 
particle correlations but also other aspects of heavy-ion 
collisions. 



B. Energy dependence of BE parameters 

The dependence of the parameters Rside, Rout, and 
Riong on transverse momentum in central heavy-ion col- 
lisions at various center of mass energies is summarized in 
Fig. 18. All the data correspond to measured pion-pion 
correlations at midrapidity, though it should be noted 
that the analyses differ in details such as the choice of 
the transverse momentum variable and the treatment of 
the Coulomb interaction. 

From Fig. 18 it is obvious that both the absolute values 
and the rate of decrease of the radii with increasing trans- 
verse momentum exhibit only a weak energy dependence. 
Since correlations are assumed to reflect the evolution of 
the system and its striicture at freeze-out, these system 
properties therefore seem to depend only weakly on the 
initial conditions of the system. As only central colli- 
sions are considered, the initial overlap region is similar, 
but the energy density increases by an order of magni- 
tude between AGS and RHIC energies. This increase of 
the initial energy density causes a pronounced change in 
other observables, e.g. a strong increase of the particle 
abundances and a change in the particle composition. 

Fig. 19 shows the energy dependence of the BE radii 
at kt — 0.2 GeV/c. The plotted values have been calcu- 
lated by linear interpolation of the data in Fig. 18. The 
parameter Rside decreases at AGS energies and remains 
almost constant at higher energies, the Rout parameter 
remains almost constant, only the Rumg parameter shows 
a trend to increase with the center of mass energy. 

In many model calculations the parameter Raids at 
kt = Q is foimd to be closely related to the transverse 
size of the system. Under this assumption, the appar- 
ent independence of Rgide on the center of mass energy 
in Fig. 19 means that the system at freeze-out reaches 
roughly the same radial extension irrespective of the ini- 
tial energy density. 

Provided Eq. 15 is a valid approximation and the 
freeze-out temperature is constant, the slight increase of 
the Riong parameter might be connected to a longer over- 
all lifetime of the pion source created in central collisions 
at larger center of mass energies. 

The ratio Rqul to Rside has been related to the emission 
duration in the emission function. According to Fig. 19 
this ratio is close to unity independent of the center of 
mass energy. This corresponds (see Eq. 17) to a small 
emission duration, i.e. to a sudden freeze-out. 



In [46] it has been suggested that the freeze-out occurs 
when the mean free path of the pions. 



reaches a value of about 1 fm. Here, N denotes the num- 
ber of particles in the freeze-out volume Vf and a the pion 
cross section. The freeze-out volume has been related to 
the radii by 

Vf = {2Trf''-Rl^,-Rl,r.g. (19) 

Fig. 19 shows this quantity in the lowest panel. In [46] a 
minimum in the energy dependence of Vf was observed 
at lower SPS energies, where the system transforms from 
a baryon- to a meson- dominated system. The radii deter- 
mined in the present analysis do not confirm this result. 

Another ansatz for the interpretation of the beam en- 
ergy dependence of the extracted radii is given in [47, 48] . 
The product of the radii Rside, Rout, and Riong is used 
to derive the phase space density of the system aver- 
aged over momentum and configuration space. Since this 
quantity is believed to be conserved at the later stages 
of the assumed isentropic evolution of the system, it is 
closely related to the conditions of the system at earlier 
times, when possibly a phase transition from deconfined 
matter to a hadronic system occurred. It was found [48] 
that the averaged phase space density shows an energy 
dependence consistent with the appearance of a decon- 
fined phase at lower SPS energies. 

C. Rapidity dependence at SPS 

This section discusses possible implications of the ra- 
pidity dependence of the BE correlation function of pion 
pairs. With the choice of a particular rapidity interval, 
particle pairs are selected that are emitted from certain 
parts of the source: Assuming an emission function as de- 
scribed in Eq. 20, high pair rapidities correspond to large 
z-values of the last interaction (i.e. emission) point. Us- 
ing relations 15 - 17, the small decrease with increasing 
rapidity of the Rside parameter at low transverse momen- 
tum implies a smaller transverse extension of the source 
as well as a higher expansion strength for the forward 
region of the system. Since Rout does not change with 
increasing rapidity, but Rside decreases slightly, the emis- 
sion duration parameter, proportional to the difference 
of the two, also increases slightly. The Riong parame- 
ter varies with rapidity, especially at low kt, but not in 
a systematic way. Hence conclusions are hard to derive 
without consulting more refined model calculations. The 
cross term Routiong clearly rises with increasing rapidity, 
indicating the strong longitudinal expansion. A different 
reference frame and parametrization as applied in [14] is 
better suited to quantify the strength of the longitudinal 
expansion. The results in fact demonstrated approximate 
boost-invariance . 
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FIG. 18: (Color online) Summary of the kt dependence of the radii at various center of mass energies at midrapidity. Data 
are from AGS [20] (squares), this NA49 analysis (circles), PHOBOS [28] (triangles at yil^ = 62 GeV), and STAR [25, 26] 
(triangles at = 130, 200 GeV). 



Experimentally a clear increase of the radii with de- 
creasing impact parameters of the collision has been ob- 
served [24, 26]. Decreasing the impact parameter leads 
to larger rapidity densities. Therefore it has been specu- 
lated that there is an explicit relation between the rapid- 
ity density and the size of the radii. Since the rapidity 
density also depends on rapidity (see Fig. 9), this quan- 
tity can also be varied - at constant beam energy - by 
changing the rapidity interval. Fig. 20 displays the de- 
pendence of the radii on the rapidity density in central 
heavy-ion collisions. Though the parameter Riong in- 
creases with increasing rapidity density, both Rside and 
Rout depend only little on rapidity density. This con- 
trasts with the impact parameter dependence, where all 
three parameters increase from peripheral to central col- 
lisions. 



D. Blast- wave model 

In this section an analytical parametrization (blast- 
wave (BW) model) by Lisa and Retiere [49] of the 
hadronic freeze-out is used to determine key characteris- 
tics of the evolution of the collision from the measured 
radii and transverse mass spectra. The parametrization 
is based on a hydrodynamical ansatz, in which the expan- 
sion of the particle source is quantified in the most simple 
case by only two parameters: The transverse expansion 
velocity and the temperature. This approach was intro- 
duced by Schnedermann et al. [50] and has been further 
developed by several groups, e.g. [51-54]. Details of the 
parametrization applied here are given in [49]. 



The model is based on the assumption of a classical 
emission function S{x,p), 

S{x,p) = mt cosh(77 - y)exp ^—^^^-^^ 

■n{r){exp{-p-u{x)/T + s)y\ (20) 

which represents the probability that a pion of momen- 
tum p decouples at the space point x from the particle 
source. The emission function in form of Eq. 20 is valid 
only in the longitudinal co-moving system. On the right 
hand side of Eq. 20 the transverse mass of the emitted 
particle is labeled mt, the space-time rapidity of the emis- 
sion point rj = ^\n[{t + z)/{t ~ z)], and the particle ra- 
pidity y. If the particle obeys Bose-Einstein statistics s is 
set to s = — 1, if Fermi-Dirac statistics has to be applied 
it is set to s = -f 1. 

The emission is assumed to occur according to a 
Gaussian distribution in longitudinal proper time t = 
\/t'^ — with the mean value labeled tq and the standard 
deviation At. The density distribution of the source is 
modeled as a uniform cylinder, i.e. ri(-\/a;2 + y^) — f2(r) 
is set to unity if r < Rgeo and zero otherwise (so called 
'box profile'). 

The Boltzmann factor exp(— p-M(a::)/T) arises from the 
assumption of local thermal equilibrium at a temperature 
T within a source element moving with four velocity u(x). 
In the transverse plane, the flow rapidity is assumed to 
increase linearly with the distance from the origin (beam 
axis), and the maximum transverse rapidity po at the 
cylinder surface Rgeo is treated as a free parameter. The 
longitudinal flow velocity is set to wl = zjt as required 
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FIG. 19: (Color online) Dependence of the BE radii, the ra- 
tio Rout/Rside and the volume Vf in central PbPb (AuAu) 
collisions on beam energy at midrapidity. The fct-dependence 
(Fig. 18) has been interpolated to get numerical values of the 
radii at kt = 0.2 GeV/c. 

by longitudinal boost-invariance. This assumption, intro- 
duced by Bjorken [55], simplifies the model, but the finite 
extension of the system and possible deviations from a 
boost-invariant expansion are neglected. 

Fixing the values of the free parameters in Eq. 20 and 
integrating it yields the predictions for the single-particle 
spectra: 

(X f d^xS{x,p). (21) 
Ptdpt J 

An optimal set of parameters can be obtained by com- 
paring these predictions with measured transverse mo- 
mentum spectra and employing a fit procedure. Here 
it should be pointed out, that the BW parametriza- 
tion neglects the fact that a large fraction of the mea- 
sured hadrons, particularly low-pt pions, stems from res- 
onance decays. The measured spectra actually repre- 
sent the sum of the spectra of the primary particles and 
the decay spectra of strongly decaying resonances. Ex- 
perimentally, these two contributions can not be distin- 
guished. Effects of resonance feeddown on fits with the 
BW model were minimized by using pion spectra only 
for mt > 0.4 GeV/c. 

In the BW model [49] the radius parameters are ob- 



FIG. 20: (Color online) Dependence of the BE radii on dn/dY 
in central PbPb (AuAu) collisions. Closed symbol correspond 
to measurements at midrapidity, squares correspond to AGS 
energies, circles to SPS energies, and triangles to RHIC ener- 
gies. Open symbols represent NA49 data at forward rapidity 
1.0 < y,r7r < 1.5. All radii have been interpolated to get 
values at kt = 0.2 GeV/c. 



tained by calculating numerically the appropriate spatial 
and temporal moments of the emission function S{x,p) 
for a given set of BW parameter values in Eq. 20. Using 
an iterative fit procedure, the optimal parameters of the 
emission function can then be determined by comparison 
of the computed radius parameters with their measured 
values and /ct-dependence. 

The fitting procedure can treat single-particle spectra 
and radius parameters simultaneously, ensuring that the 
finally extracted emission function is better constrained 
than in studies which consider only one set of observables. 

In the present application of the model an azimuthally 
symmetric source is assumed, since only central collisions 
are studied. With this restriction the model features five 
free parameters: the freeze-out temperature T, the max- 
imum transverse rapidity poj the freeze-out transverse 
radius Rgeo, the overall lifetime of the source tq and the 
emission duration At. Elliptic flow measurements are 
not included in this study, since the analysis focuses on 
central collisions. In the following, this parametrization 
will be fitted to radius parameters and to transverse mass 
spectra at midrapidity using only statistical errors. 

The data points in Fig. 21 show single-particle spectra 
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measured by NA49 at five beam energies at midrapidity 
and the lines indicate the fit results. The rrit dependence 
of the production of the various particle species is well 
described by the parametrization. Since the model only 
aims to reproduce the shape of the spectra, the normal- 
isation, i.e. the particle yield, was treated as a free pa- 
rameter in the fitting procedure. Solid lines in Fig. 21 are 
obtained from the fitting procedure when negative pions, 
protons and the radii are used as input. The predictions 
for K"*", K~ and antiproton spectra were then calculated 
using the fitted parameters and are shown by the dashed 
curves in Fig. 21. They are also in excellent agreement 
with the measurements. 

Fig. 22 shows the kt dependence of the parameters 
Rside, Rout, and Riong at different incident beam ener- 
gies at midrapidity compared with the fit values. Both, 
the absolute values and the transverse momentum de- 
pendences are well reproduced. Only the gradient of the 
fct-dependence of the parameter Rside tends to be smaller 
in the parametrization than in the measurement. 

The energy dependence of the extracted BW model pa- 
rameters is plotted in Fig. 23. In addition to results ob- 
tained from fits to NA49 data at SPS energies, the figure 
also includes parameters derived by the same fit proce- 
dure from data at lower energies at the AGS [20, 56, 57] 
and at higher energies at RHIC [26, 58-60]. From ex- 
periments at the AGS radii for pions and transverse mo- 
mentum spectra of only pions and protons are available. 
For consistency, also the fits at higher energies used as 
input only radii and transverse momentum spectra of pi- 
ons and protons. Open symbols in Fig. 23 represent the 
result when only pion and proton spectra are fitted and 
two-particle correlations are disregarded. In this special 
case only two parameters are determined, T and po- 

The temperature parameter tends to increase with the 
available center of mass energy, whereas the surface ex- 
pansion rapidity po saturates at lower SPS energies. The 
kinetic freeze-out temperature of about 80-100 MeV is 
lower than the hadro-chemical freeze-out temperature at 
an earlier stage of the collision. The latter can be de- 
rived from the particle species composition resulting in 
values of about 140-160 McV [61]. This observation com- 
plies with the assumption that the system cools down by 
rapid expansion first through chemical freeze-out (when 
inelastic collisions stop in the fireball) and then through 
kinetic freeze-out (when also elastic interactions cease). 

The transverse radius at freeze-out is independent of 
the beam energy at about 12 fm, which is approximately 
twice the size of a lead nucleus, indicating the growth of 
the system. The lifetime of the source depends weakly 
on energy; it is about 4-6 fm/c at lower energies and 
increases to 8 fm/c at RHIC energies. 

The emission duration takes a finite value at lower en- 
ergies and increases up to the top SPS energy, then it 
saturates at a value of ~ 3 fm/c. 
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FIG. 21: (Color online) Transverse mass spectra of pions, 
kaons, and protons in central PbPb collisions at midrapidity 
measured by NA49. mo is the particle rest mass. Lines indi- 
cate the fit results, numerical values of the fit parameters T 
and po are also shown. Further details are given in the text. 
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FIG. 22: (Color online) fct-dependence of the radius param- 
eters in central PbPb collisions at midrapidity (data points) 
compared to fit results (lines). 



E. Comments on other models 

The Buda-Lund collaboration uses a parametrization 
similar to the BW model described in the previous sec- 
tion to describe simultaneously single-particle spectra 
and radii. The parametrization is introduced in detail 
in [51, 62] and has been successfully appHed to SPS and 
RHIC data [63, 64]. The main extensions compared to 
the simpler blast-wave model are the assumption of a lo- 
cal temperature T{x) instead of a single global tempera- 
ture T, and the introduction of a local chemical potential 
fj-{x). Preliminary results using this parametrization are 
similar to those obtained with the BW model discussed 
in the previous section. A detailed analysis of the data 
presented here is in progress. 

A different approach is provided by pure hydrodynamic 
calculations, an idealized but well-defined limiting case 
scenario for the evolution of the system generated in cen- 
tral heavy- ion collisions. This kind of model relies only 
on the assumption of an equation of state, a set of initial 
conditions, and on a specific freeze-out criterion. The 
matter is treated as an ideal, locally thermalized fluid 
whose dynamics are governed by the relativistic hydro- 



dynamic equations. The freeze-out is usually modeled to 
occur at a fixed temperature. If the local temperature 
of a fluid cell drops below the freeze-out temperature, 
the properties of this cell are converted into an ideal gas 
of hadrons and hadronic resonances (Cooper-Frye pre- 
scription) [65, 66]. This leads to an emission function 
which can be used to calculate the two-particle correla- 
tions. This class of models successfully describes bulk 
properties like the elliptic flow and the transverse mo- 
mentum spectra of pions. However, the models often fail 
to reproduce the experimentally observed fc(-dependence 
of the radii [67, 68]. The overall lifetime in this class 
of models is usually found to be large, hence the Riong 
parameter is often overpredicted. Depending on the as- 
sumed equation of state and initial conditions, the Rside 
parameter is often underpredicted and Rout is overpre- 
dicted. Consequently, the ratio of these two parameters, 
related to the emission duration, is not well reproduced. 
If the system passes through a first-order phase transi- 
tion the Rout parameter was predicted to become much 
larger than Rside [69] . Experimentally this signature has 
not been observed, as can be seen from Fig. 19. How- 
ever, more refined models are able to describe the data. 
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FIG. 23: (Color online) Energy dependence of the model pa- 
rameters T, po, TO, At, and Rgeo- Input to the fit procedure 
are the spectra of n~ and protons, and the radii. Open sym- 
bols correspond to a fit of only the spectra, disregarding the 
radii. Further details and references are given in the text. 
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E.g., in [30] a protracted surface emission is introduced 
instead of the conventional Cooper- Frye prescription [66] . 
In this case the fcj-dcpendcncc of the ratio Rout / Rside is 
well reproduced. The failure of pure hydrodynamic mod- 
els might also be due to the fact that the correlations are 
determined at the final stages of the evolution, when the 
mean free path of the particles is finite, possibly of the 
order of the size of the system. In contrast, the basic as- 
sumption of the hydrodynamic approach is a zero mean 
free path. 

A different concept is pursued by microscopic models. 

These provide an approach complementary to the hy- 
drodynamics inspired ansatz described above. In trans- 
port approaches, the heavy-ion collisions are modeled by 
initial particle production from string fragmentation fol- 
lowed by evolving these particles along straight-line tra- 
jectories punctuated by collisions according to free-space 
cross-sections. The result of this microscopic simulation 
is the state of the system at freeze-out, i.e. the momenta 
and space-time points of the last interaction. A well 
known implementation of this kind of models is UrQMD, 
which is described in detail in [70-72]. Quantum statis- 
tics particle correlations are in principle not included in 
transport models. Therefore, BE correlations have to be 
introduced in a second step, e.g. by using an afterburner 
hke the CRAB algorithm [73, 74]. 

The UrQMD model combined with the CRAB algo- 
rithm was compared recently to the radius parameters 
measured at SPS energies [75]. The model is able to 
reproduce qualitatively the kt and rapidity dependence 
of the radii at all energies. In particlular it reproduces 
the weak rapidity dependence of the parameters Rside, 
Rout, and Riong, and the strong increase of the cross term 
Routiong with increasing forward rapidity. The parame- 
ter Rside is slightly underestimated and Rout is overesti- 
mated at larger transverse momenta. Consequently, the 
ratio Rout/ Rside exceeds the measured values. 

VIII. CONCLUSION 

This paper presents a detailed study of 7r~7r~ Bose- 
Einstein correlations in central heavy-ion collisions at five 
SPS energies. The kt and rapidity dependence of the 
fit parameters A, Rside, Rout, Riong, and Routiong is ob- 
tained. A decrease of the radius parameters with increas- 
ing transverse momentum is observed. This behaviour 



has been reported before and is usually attributed to the 
rapid expansion of the pion source. 

Absolute values as well as the gradient of the trans- 
verse momentum dependence do not change with incident 
beam energy. A blast-wave parametrization of the evolu- 
tion of the system motivated by hydrodynamics yields a 
freeze-out temperature of 80-100 MeV, a transverse ex- 
pansion velocity of the order of 72 % of the speed of light 
at the surface, a transverse radius of the system of about 
12 fm, a total lifetime of about 6 fm/c and a finite emis- 
sion duration of about 3 fm/c. The energy dependence of 
the radius parameters as well as of the fitted blast-wave 
model parameters shows no particular structure. 

The study of the rapidity dependence of BE correla- 
tions reveals that only the Rside parameter depends sig- 
nificantly on the selected pair rapidity. This parameter 
decreases systematically from midrapidity to forward ra- 
pidity. Rout depends only weakly on rapidity and Riong 
shows some variations at low transverse momenta, but 
not in a systematic way. The weak rapidity dependence 
of the radii contrasts with the strong change of the ra- 
pidity density of produced particles. 

The data set presented here provides stringent con- 
straints for models trying to describe the evolution of 
central heavy-ion collisions at SPS energies. Experiment 
NA61, an upgrade and continuation of the NA49 experi- 
ment, will extend the systematics with measurements of 
smaller collision systems [76]. Moreover, an analysis of 
BE correlations in NA49 using the source imaging tech- 
nique [39] is in progress. 
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